Detecting relics of a thermal gravitational wave background in the early Universe 
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Abstract 

A thermal gravitational wave background can be produced in the early Universe if a radiation dominated epoch precedes 
the usual inflationary stage. This background provides a unique way to study the initial state of the Universe. We discuss 
the imprint of this thermal spectra of gravitons on the cosmic microwave background (CMB) power spectra, and its possible 
detection by CMB observations. Assuming the inflationary stage is a pure de Sitter expansion we find that, if the number 
of e-folds of inflation is smaller than 65, the signal of this thermal spectrum can be detected by the observations of Planck 
and PolarBear experiments, or the planned EPIC experiments. This bound can be even looser if inflation-like stage is the 
sub-exponential. 
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I. INTRODUCTION 



Understanding the expansion history of the Universe is a fundamental task of modern cosmology. The current 
observations of the cosmic microwave background radiation (CMB) pj, large scale structure 0], Type la supernova 
amongst others, have provided us with a relatively clear picture of the expansion history of the Universe since 
the photon decoupling at the redshift z ~ 10 3 . The history of the Universe prior to decoupling is currently primarily 
deduced through indirect evidences, such as the primordial abundances of the light elements. The primordial abun- 
dances hold the information about expansion history of the Universe at the epoch of nucleosynthesis corresponding 
to redshift z ~ 10 10 . 

Going back even further, in order to solve the flatness, horizon and monopole problems in the standard hot big-bang 
cosmological model, various inflation- like scenarios [4j (or the so-called 'kick' scenarios in [5() for the expansionary 
history of the Universe at very high redshift have been proposed. These scenarios generically predict a nearly scale- 
invariance primordial power spectra of density perturbations and gravitational waves [||, and have been indirectly 
supported by the observation of the large scale structure and CMB temperature and polarization anisotropies power 
spectra [H, 0|. 

In this Letter we shall analyze the possibility of observing imprints from a pre-inflationary stage. These imprints 
could give us a glimpse at the physical conditions of the very early Universe right to the time of its birth. It is natural to 
suppose that the Universe underwent a radiation dominated stage of expansion prior to the inflation-like acceleration 
phase [f\. Furthermore, it is reasonable to assume that, during this stage at temperatures higher than ~ 10 19 GeV 
a thermal equilibrium between the various components, including gravitons, is maintained through gravitational 
interaction. In this scenario, as the Universe cooled down and the gravitons decoupled, a background of thermal relic 
gravitons with a black-body spectrum would be left behind 0, Q . This thermal background of gravitational waves 
will garner uncontaminated information about the thermal pre-inflationary period, and it's detection would give a 
unique chance to probe the physics of pre inflationary Universe inaccessible by other means. 

As the Universe expands, and undergoes inflationary expansion, the gravitational waves would be strongly redshiftcd 
to very low frequencies fioj . The thermal peak frequencies would correspond to the range probed by CMB experiments 
v oj 10 -18 — 10~ 15 Hz. The gravitational wave background could leave an observable imprint in the temperature and 
polarization anisotropies of the CMB, which is expected to be detected by observations in the near future. In the 
current work, we shall discuss this signature in the CMB, and analyze the possibility of detecting this signature by 
the upcoming observations. 

In Section [TT] we shall briefly consider the main physical motivations for the existence of thermal gravitational 
wave background. We give simple estimates for the main characteristics of this field and relate them to inflationary 
parameters. In Section IIHI we analyze the affect of this background on the CMB temperature and polarization 
anisotropies. Based on the WMAP 5 year data, we place upper bounds on conformal temperature (explained below). 
We compare the expected signal with the sensitivity of the upcoming CMB experiments, and study the feasibility 
of detecting this background. Finally, in Section IIV1 we discuss the physical implications of an observable thermal 
gravitational wave background. 

Throughout this Letter, we will work with units in which c = h = ks = 1. 

II. THERMAL BACKGROUND OF RELIC GRAVITATIONAL WAVES 

Let us discuss the possibility of generating black-body spectra of gravitational waves at very high energy scales in 
the early Universe. We shall assume that the Universe was radiation dominated before the inflationary epoch, and 
that all the particle species were highly relativistic. The interaction rate for particles interacting solely through the 
gravitational force would be T ~ T 5 /Mi, where T is the physical temperature in the Universe, and M p \ = G -1 / 2 = 
1.22 x 10 19 GeV is the Planck energy [1, [ll|. So long as this interaction rate is large in comparison with the expansion 
rate characterized by Hubble parameter H ~ T^/M p i, the gravitons would remain in thermal equilibrium with other 
particles. However, this equilibrium will be violated once T < H . At this time, the gravitons would decouple from the 
other particle species, leaving behind a free-streaming thermal graviton background. At the time when the gravitons 
decoupled, the interaction rate is T ~ H, i.e. T ~ M p i. It should be mentioned that the conditions for a thermal 
equilibrium of gravitons with other particles existed only at very early times T > M p i, where the physics is not 
properly understood. For this reason, the graviton background might not have a black-body spectrum, although it 
would be in a non- vacuum state. However, in this Letter, we shall restrict our considerations to a black-body spectrum 
as a concrete example of a non- vacuum state, similar to considerations in [1, Q. We shall assume that this thermal 
graviton background decoupled at temperature T ~ M p i from the surrounding matter. 

With the expansion of the Universe this graviton background would preserve its thermal spectrum, but would be 
strongly rcdshifted to very low temperatures. It is instructive to estimate the temperature of this background at the 
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FIG. 1: The evolution of Hubble radius, H 1 , and various physical sizes with the expansion of the Universe. The left panel 
shows the result in the scenario with Ah > A p ;, and right panel shows the result in the scenario with Ah < A p ;. 



present epoch. In order to proceed, it is convenient to separate the history of the Universe into three stages: the 
initial radiation dominated stage, the inflationary stage and the post-inflationary stage. In terms of this division, the 
present day temperature 7o of the graviton background is 

Mp\ CLi n f a p do 

where a, , cij n / , a p and ao are the values of the scale factor at the time of graviton decoupling, the beginning of the 
inflationary stage, the end of inflation (beginning of the post-inflationary stage) and the present day, respectively. 
During and after the inflationary stage, the temperatures of the graviton field and that of the rest of the particle 
species behave in significantly different manner. The temperature on the graviton field strictly decreases with the 
expansion. On the other hand, at the end of inflation the temperature of the thermal bath, containing the rest of 
the particle species is significantly boosted by the process of reheating. Assuming that the observed CMB is the relic 
of this thermal bath, its temperature T2 at the beginning of the post-inflationary stage can be related to value of 
the scale factor at this stage through the relation a p /a^ = (2. 73K/72)(3. 91/106. 75) 1//3 [g|. Using this expression, and 
denoting the temperature of the thermal bath at the beginning of the inflationary stage as 7{ , the equation (fT]) can 
be rewritten as 

T ~ 8.0 x 10- 27 (T 1 /T 2 )e 60 - JV K, (2) 

where TV = log(a p /aj n /) is number of e-folds during the inflationary stage. This spectrum is peaked at the frequency 

i/~4.7x 10- 16 (f 1 /f 2 )e 6 °- A, Hz. (3) 

As can be seen, the peak frequency depends on the value of the number of e-folds during inflation, and the ratio of 
the temperatures T\ and 7^. 

The number of e-folds N cannot be very small. It must be larger than some iV m i n to account for the isotropy 
and homogeneity of the observed Universe. If we take the inflationary stage to be an exact de Sitter expansion at 
7~i ~ 10 16 GeV, a value N m i n ~ 60 follows (see for instance [HI, [l3j]). In some realistic inflationary models, this 
minimum value could be even lower N m [ n ~ 46 [l2T| , On the other hand, it is also important to consider the upper 
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limits on the e-folds number. In general, the value of N does not have a necessary upper bound. However, if we focus 
on specific inflationary models, the value of iV cannot be too large. For example, in order to account for the observed 
scalar spectral index n s = 0.96, in the inflationary model with the potential form V(4>) — A 4 (</>//i) p with p — 2, one 
has N — 50. And for the model with p = 4, one has N — 74 12]. In this section, as a rough estimation, we shall 
restrict our analysis to a typical range for the number e-folds N £ (60, 70). 

In the case of an inflationary stage characterized by de Sitter expansion one has T\ ~ T^,. Assuming N £ (60, 70), we 
get v £ (2.1 x 10 -20 , 4.7 x 10 -16 ) Hz. In the case of an inflationary stage with sub-de Sitter expansion, the value of 71 is 
typically larger than that of T 2 . Assuming f x /% = 100 and N £ (60, 70), we obtain v £ (2.1 x 10~ 18 , 4.7 x 10~ 14 ) Hz. 
It is important to note that in both the cases the peak frequency is in the range typically probed by the observations 
of the CMB. It is therefore reasonable to look for the signature of the relic thermal graviton spectrum in temperature 
and polarization anisotropies of the CMB. In the following section we shall address this question in more detail. 
However, it is also worth noting that, when the e-folds number is very large, N ~ 100 for example, the peak frequency 
will be v ~ 10 _30 Hz, leaving no hope for observing the thermal nature of the graviton background. 

A natural question arises in the context of the evaluated above typical peak frequencies of the thermal graviton 
spectrum. Namely, is it possible for the peak wavelength X p i of the thermal graviton spectrum to be larger than 
the minimum scale Aj n / at which the Universe is homogenous and isotropic? The minimal scale of homogeneity and 
isotropy is set by the Grishchuk-Zel'dovich effect [3] at Am/ ~ 500 Xh, where Ah = Hq 1 is the present day Hubble 
radius. If it were the case that A p ; > Xinf, the thermal nature of the graviton background would not be observable, 
since it would correspond to a peak frequency v < 10 _21 Hz, smaller than those accessible to CMB measurements. 
However, as illustrated in Fig. [T] such a situation does not arise. In Fig. [TJ we plot the evolution of the Hubble radius 
H^ 1 (black solid line) through the different stages of expansion of the Universe. During the radiation-dominated 
stage (before and after the inflationary stage), H' 1 cx a 2 , and in the following matter-dominated era H^ 1 cx a 3 / 2 . 
In the inflationary stage, ff -1 is nearly a constant. The physical length associated with a length-scale of a fixed 
conformal length scales as L cx a with the expansion. In Fig. [T]we show three such length scales: A,,; is the Planck 
length scale corresponding to the peak in the thermal graviton spectrum, Xh is the present day Hubble radius, and 
Xi n f is the scale associated with large scale homogeneity and isotropy. Fig. [1] shows that, if Aj n / is taken to be 
equal to the Hubble length at the end of initial radiation dominated era at the beginning of the inflationary epoch, 
then it necessarily follows that Aj n / > X p i. This is a general statement, relying only on the assumption that the 
pre-inflationary state satisfies an equation of state p/p > —2/3. Thus, there are two possible scenarios. In the first 
case, when Xh > X p i (left panel in Fig. [1]), corresponding to a small number of e-folds, the thermal nature of the 
graviton spectrum would be observable with CMB measurements. In the second case, Aj n / > X p i > Xh (right panel 
in Fig. [1]), corresponding to a large number of e-folds, the thermal nature of the graviton spectrum would not be 
observable with CMB measurements. 

Before proceeding to the analysis of the observable signature of a thermal background of gravitational waves in 
the anisotropies of the CMB, let us briefly analyze the resultant primordial power spectrum of this backround. The 
term primordial power spectrum denotes the power spectrum of relic gravitational waves in the radiation dominated 
epoch after the end of inflationary epoch when the wavelength of interest are significantly larger than the horizon. 
The thermal nature of the gravitational wave background prior to the inflationary stage will leave a distinct signature 
in the primordial power spectrum, which will then translate into a specific features in the power spectrum of CMB 
anisotropies. We start with the gravitational field of a slightly perturbed Friedmann-Lemaitre- Robertson- Walker 
universe given by 



where r) is conformal time, Sij is the Kronecker delta symbol, and the metric perturbation field hij only contains 
contribution from pure gravitational waves. The gravitational wave field has two modes of polarization, which can be 
expanded over spatial Fourier harmonics 



ds 2 = a 2 (ri){-dri 2 + (<% + h lj )dx i dx j ], 



(4) 




ikx 



8=1.2 




(5) 



where pff (k) is the polarization tensor [l5[ 



. The power spectrum is defined as 



(* k /£,> = ^P t (*)« 3 (k-k'), 



(6) 
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where the angle brackets indicate an ensemble average (see for example [la]). Assuming vacuum initial conditions at 
a pre-inflationary stage, the Fourier coefficients satisfy the relations 

(4c k ,)=<5 3 (k-k'), (c k c k ,) = <44')=0. (7) 

However, if we assume that the gravitational wave field was in thermal equilibrium, then the first relation in ([7]) 
modifies to [j| 

<4^>= (l + ^ TT )'5 3 (k-k'), (8) 

where T is the conformal temperature of the gravitational wave background. The conformal temperature is related 
with the present day physical temperature Tq by the relation 

T = %a . (9) 

In the present Letter we set a = 1, so that T = T . 

The evolution of the universe through the inflationary expansion phase converts the initial vacuum into a multi- 
particle squeezed vacuum state characterized by a power law primordial spectrum [l6j |. In the case of an initially 
thermal background, the primordial spectrum takes a modified form [TJ, [TtJ : 



P t {k) = A t (k ) 







' k ' 


Go, 


1 coth 






2T 



(10) 



where n t is the spectral index, which is close to zero for typical inflationary scenarios. A t (ko) is the amplitude of 
the spectrum at the pivot wavenumbcr fco- For large wavenumbers k 3> T, the power spectrum Pt(k) oc k' H is 
indistinguishable from the case of vacuum initial condition. However, for small wavenumbers k <C T, the spectrum 
shows dissimilarity Pt(k) oc k nt ~ l . Thus, generically, a thermal background exhibits difference in the power spectrum 
to the initial vacuum background at low wavenumbers k <T . This point was previously raised in [l8l |. were the author 
argued that modifications of Pt, due to an initial thermal spectrum, could be constrained by CMB observations. In 
the following discussion, we will also discuss the constraint on the conformal temperature T by the 5-year WMAP 
observations. 



III. THE IMPRINT IN THE CMB AND ITS DETERMINATION 



In this section we shall analyze the observable signatures of the thermal background of relic gravitational waves 
in the power spectrum of temperature and polarization anisotropies of the CMB. The power spectrum of CMB 
anisotropics due to gravitational waves are evaluated by solving the radiative transfer equation in the framework 
of perturbation theory [l9|, [2(J [2l| (see also [HI). In general, gravitational waves leave their imprints in the four 
anisotropy spectra Cf x '(XX' = TT, TE, EE, BB) HUES]. However, in this work, we shall restrict our discussion 
to the power spectrum of the so-called i?-mode of polarization Cf B . The B-mode is solely generated (neglecting 
the possible foregrounds) by gravitational waves, and thus provides a clean channel for their detection. The features 
in the gravitational wave power spectrum P t at wavenumber k translate predominantly into features in the CMB 
power spectrum at multipole i ~ k x 10 4 Mpc [23j]. For this reason, we can hope to detect the signature of thermal 
gravitational wave background in the CMB for temperatures T > O.OOOlMpc -1 . In Fig. [21 we plot the CMB power 
spectrum Cf B for various values of T. As expected, the signature of the thermal background is predominantly located 
at multipoles I < T x 10 4 Mpc. 

The current CMB observations have yet to detected the B-mode of polarization. The 5-year WMAP observation 
give an upper limit for the -B-mode at lower multipoles, with the average value £(£+ l)Cf B /2ir at I = 2, 3, 4, 5, 6 being 
smaller than 0.15/xK 2 (95%C.L.)[24|. This upper limit, allows to place constraints on on the value of the conformal 
temperature depending on the tensor-to-scalar ratio r, shown in Fig. [3] The tensor-to-scalar ratio r characterizes 
the overall contribution of gravitational waves to the CMB anisotropies, and is defined as the ratio of the primordial 
power spectra of gravitational waves to density perturbations (see for example [23jj ) . As expected, the constraints on 
T scale inversely proportional to the value of r. For r = 0.1 the constraints read T < 0.016Mpc _1 , while for r = 0.3 
we have T < 0.005Mpc _1 . 

Let us now turn our attention to the analysis of prospects for future CMB observations. In order to determine 
the potential of the future CMB observations to constrain the conformal temperature of the thermal gravitational 
wave background we shall use an approach based on the Fisher information matrix. In terms of the Fisher matrix, 
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FIG. 2: The B-mode polarization power spectrum Cf B in the models with different conformal temperature T. In all the 
spectra nt — and r = 0.1 at the pivot wavenumber ko — 0.2Mpc -1 . The cosmological parameters are set at the WMAP5 
best-fit value [J]. 

the precision on a parameter pi that can potentially be attained is given by A Pl ~ (F^ 1 )*/ 2 [HHil. In the present 
analysis we shall restrict to two free parameters T and r. 

The other cosmological parameters can be determined with high precision from the other CMB spectra Cj T , Cf E 
and Cf E , and their inclusion in the analysis will not significantly alter the results. Therefore, we shall fix the other 
cosmological parameters at their fiducial values fib = 0.0456, f2 c = 0.228, J1a = 0.726, Qk = 0, h = 0.705 and 
n s — 0.96, A s — 2.036 x 10~ 9 at the pivot wavenumber fco = 0.2Mpc _1 For simplicity, in this Letter, we shall also 
fix the gravitational wave spectral index rit = 0, corresponding to a flat (scale-invariant) power spectrum. 

The Fisher matrix can be written as I2a| 



Pir BB ar BB 1 

*«-2^ dPi dPj (AD?*r 



where i,j — 1,2 correspond to the parameters T and r respectively. The quantity ADf B is the standard deviation 
of the estimator D? B O 



AD ™ = iWT^ {cr + Nr) - (12) 

In the above expression the noise power spectrum N BB and the sky cut factor / s k y are determined by the specificities 
of the particular CMB experiment, see for example [la |. 

We now estimate the potential to constrain the conformal temperature T based on combining the data from the 
space-based Planck satellite [27j and the ground-based PolarBear experiment (28[. The former is sensitive to the 
B-mode of polarization at I < 20, while the latter one is sensitive at I > 20. Their combination provides an excellent 
opportunity to determine the gravitational wave signal [29l | . The corresponding instrumental noises and the sky cut 
factors can be found in [2?], [Hj]. In addition to the instrumental noises, in calculating the noise power spectrum N BB , 
we have included the contribution due to the cosmic lensing effect [30(. 
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tensor-to-scalar ratio: r 



FIG. 3: The upper limit of the conformal temperature T as a function of the tensor-to-scalar ratio r, based on the WMAP5 
observations. 

Using expression we calculate the value of AT for given values of parameters T and r. For a fixed value of 
r, the smallest measurable value of T is determined by the condition T = AT. This lowest bound on the detectable 
signal is shown by blue line on Fig. 2] For a a typical value r — 0.1, the attainable limit on the conformal temperature 
is T = 1.8 x 10 _4 Mpc _1 . The lower limit on the detectable signal could be further improved by the proposed EPIC-2m 
experiment [3l[ shown with the red line on Fig. @] In this case, the lower limit is T = 1.4 x 10~ 4 Mpc -1 for values 
r > 0.005. 



IV. DISCUSSION AND CONCLUSIONS 

The possible determination or placement of upper limits on the conformal temperature of the thermal gravitational 
wave background would allow to place interesting constraints on the physics of the inflationary era. To see this, let 
us rewrite Eq. ^ using Eq. ^ in the form 

T = 0.017 x e 60 ^ Ar (t 1 /t 2 )Mpc- 1 . (13) 

Assuming a fiducial value r = 0.1, the 5-year WMAP data places an upper limit on the conformal temperature 
T < 0.016 (see Fig. [3]). Using Eq. (|13p. this limit translates into a constraint on the parameters of the inflationary 
expansion N > 60 + log^/Tjj). Assuming Ti ~ T2, we obtain a constraint on the number of e- folds N > 60. It is 
worth pointing out that, this bound is consistent with the e-folds parameter required to solve the flatness, horizon 
and monopole problems in the standard hot big-bang cosmological model [131 ] . 

The EPIC-2m experiment would be able to determine the conformal temperature for a broad range of tensor-to- 
scalar ratios r > 0.005 down to a limit T > 1.4 x 10 Mpc -1 (see Fig. 2]). A positive detection of non-zero conformal 
temperature by the EPIC-2m experiment would place upper bounds on the number of e-folds N < 65 + log(7i/72), 
which would depend on the ratio of background temperatures before and after inflation T1/T2. For T\ ~ T2 the upper 
bound would be N < 65, while for T\ ~ 100T 2 the upper bound would become N < 69. 
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tensor-to-scalar ratio: r 



FIG. 4: The curves show the smallest value of conformal temperature T which can be determined by CMB experiments. The 
blue line shows the sensitivity of Planck+Polarbear experiment, and the red line shows the sensitivity of the proposed EPIC-2m 
experiment. 

The above considerations show that an observable thermal gravitational wave background may allow to place 
stringent constraints on the range of viable inflationary models. More over, its detection may shed light onto quantum 
gravity effects, which become important at Planck energy scales. On the other hand, an absence of observational 
indications of a thermal background would indicate one of the two possibilities. Either the initial state of the 
gravitational wave background was not thermal, or alternatively, that the number of e-folds was N > 69 so that the 
present day conformal temperature is redshifted to T < 1.4 x 10 _4 Mpc _1 . 

Finally, it is worth pointing out that, along wi th g ravitational waves, a thermal spectrum of density perturbations 
may have also existed in the very early Universe [32l. l33l. |34| . However, the nature of this spectrum depends crucially 
on the content and state of the matter in the very early Universe, questions which are yet to be fully understood. In 
addition, the evolution of the spectrum depends strongly on the physical specificities of inflation, which are also not 
fully understood. For this reason, without further simplifying assumptions, density perturbations cannot be directly 
used to probe the initial state of the Universe. 
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